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Large scale features of a randomly isotropically forced incompressible and unbounded rotating 
fluid are examined in perturbation theory. At first order in both the random force amplitude and 
the angular velocity we find two types of modifications to the fluid equation of motion. The first 
' correction transforms the molecular shear viscosity into a (rotation independent) effective viscosity. 

The second perturbative correction leads to a new large scale non-dissipative force proportional to 
the fluid angular velocity in the slow rotation regime. This effective force does no net work and 
alters the dispersion relation of inertial waves propagating in the fluid. Both dynamically generated 
corrections can be identified with certain components of the most general axisymmetric "viscosity 
tensor" for a Newtonian fluid. 
PACS: 47.27-i, 47.32.-y 
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>> | I. INTRODUCTION 

I 

The special features of turbulence in the presence of rotation have attracted the interest of many authors [1-4]. 
Relying on some experiments, the methods of study used have ranged from analytic approaches to numerical simula- 
r/3 . tions [2,5-7]. The central theme in rotating fluids is the effect of the Coriolis force, which induces anisotropy (there is 
a preferred direction, that of the rotation axis). This anisotropy is extreme in the limit of fast rotation, which actu- 
ally forces the flow to become two-dimensional (Proudman- Taylor theorem) [1]. In this work we apply perturbation 
theory to the randomly forced Navier-Stokes equation with Coriolis force as a model for the turbulent regime of a 
rotating fluid. The perturbative study of the ordinary randomly forced Navier-Stokes equation, in combination with 
the renormalization group (as an improvement of perturbation theory), has a long tradition [8-10]. 

The addition of the Coriolis force, induced by the rotation of the fluid, introduces one additional parameter, the 
angular velocity Q, or, in dimcnsionless form, either the Rossby or Ekman numbers (in addition to the Reynolds 
number) [1]. Let us focus on the Ekman number, Ek = v/(Q,L 2 ), that depends on the viscosity v and a scale L, 
roughly associated with the size of the fluid system. The Ekman number gives the relative importance of the viscosity 
and Coriolis forces. We will assume henceforth that for small il (< v/L 2 ) the turbulence is isotropic and the only 
relevant parameter is the viscosity. In this limit, the results of the study of the ordinary randomly forced Navier-Stokes 
C^) • equation hold (the random force is always assumed isotropic). 

For larger SI we will encounter new features. In fact, the only restriction on perturbative correction terms is that 
they respect the basic symmetry of the equations, in our case, the axial symmetry about the rotation axis. We will 
see that perturbation theory generates new terms fulfilling these symmetry constraints. Therefore, one must find the 
complete set of allowed terms that can arise in perturbation theory. We will determine all the terms that can be 
represented by the components of an axisymmetric "viscosity tensor" and, in particular, the ones that arise at first 
^~>j order in perturbation theory. The part of this "viscosity tensor" that is pair-antisymmetric in the indices plays a 
i-G . significant role; however, it does not lead to dissipation, and therefore, is not truly viscous. 

As in homogeneous and isotropic turbulence, we assume that the physical region of study is sufficiently far from the 
surfaces, where the boundary conditions are imposed, for them not to have any direct effect, except the presence of 
the scale L. In contrast to ordinary turbulence, this condition only implies that we can have homogeneous turbulence 
but, due to rotation, it cannot be isotropic. It is pertinent to mention here that the possibility of anisotropic forced 
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turbulence, and precisely with axial symmetry, has already been considered [11]. In this reference, however, the 
authors assume that the breakdown of isotropy occurs through a random force whose two-point correlation function 
depends on the anisotropy vector n. They derive a renormalized force proportional to second and fourth powers of n. 
In our case, we will see that the first perturbative correction is linear in SI, like the Coriolis force itself. 

This paper is organized as follows. In Section II we introduce the randomly forced hydrodynamical equations with 
rotation. We assume that the fluid is incompressible and show how to formulate them as a problem of homogeneous 
but anisotropic incompressible turbulence. Fourier analysis of the turbulent velocity field is used to organize the 
perturbation expansion [10] in Subsection II B. We introduce in Subsection II C the linear response function. Unlike the 
isotropic case, the Coriolis term leads to a non-symmetric linear response function matrix. In Subsection II D we define 
the non-linear response function and present its perturbative expansion (slow rotation). We also compute the first 
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order perturbative correction to the response function, which allows the identification of the (rotation independent) 
effective shear viscosity (proportional to the cube of the Reynolds number), and a new anisotropic force. In Section III 
we write down the most general axially symmetric "viscosity tensor" , as the existence of a preferred direction, singled 
out by the fluid rotation, requires the introduction and use of axisymmetric tensors. This rank four "viscosity 
tensor" expresses the proportionality between the fluid stress tensor and the rate of strain tensor. In isotropic and 
homogeneous incompressible turbulence, the viscosity tensor depends only on one parameter, the fluid shear viscosity. 
In the case of rotating turbulence, and for slow rotation, we find that the axisymmetric "viscosity tensor" depends 
on two parameters: the molecular shear viscosity (coming from the isotropic terms of the "viscosity tensor") and a 
new one, that arises from the anisotropic terms in the "viscosity tensor" . We also show that this new parameter 
can be identified as the coefficient of the anisotropic force calculated perturbatively (in the previous Section II D). 
Having thus established the equivalence between the perturbatively corrected randomly forced Navier-Stokes equation 
with Coriolis force on the one hand, and a (Newtonian) rotating incompressible fluid with an effective axisymmetric 
"viscosity tensor" on the other, we proceed, in Section IV, to discuss some physical consequences of the new terms 
in the perturbed fluid equations. In Subsection IV A we consider the quasi-local force induced by the anisotropic 
components of the "viscosity tensor" and show that it is proportional to the cube of the Reynolds number, and that 
it does not lead to dissipation. In Subsection IV B we study the dynamical effects of this force on the propagation of 
inertial waves. We end by discussing our results and proposing further work on the problem of rotating turbulence. 
In Appendix A we introduce the diagrammatic representation of the exact Navier-Stokes equation and the diagram 
encoding the first order correction to the response function and in Appendix B we present the technical details needed 
to carry out the perturbative calculation in the slow rotation limit. 



II. BASIC EQUATIONS AND PERTURBATION THEORY 
A. Equations of motion with random force in a rotating frame 



We start from the hydrodynamical equations for a fluid with density field p(x, t), velocity field u(x, t), pressure 
p(x,t) 7 and molecular shear (tangential) and bulk kinematic viscosities v and k, respectively. We assume that the 
fluid is rotating with constant angular velocity fi along the z axis and that it is subject to an isotropic random forcing 
per unit mass f . The mass and momentum conservation equations are 



|+V. ( pu) = 0, 

+ (u • V)u = — Vp + iA7 2 u 

at p 
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V(V ■u)-2fiAu-OA(fiAx)+f, 



(la) 
(lb) 



where d is the number of space dimensions. The dimension of space d will be kept as a free variable, although when 
we consider rotation-dependent expressions, these must be evaluated for d = 3. The momentum equation (lb) is 
supplemented with a random stirring force that leads to a statistical distribution for the velocity field and can be used 
to model turbulent flows just as is done for isotropic randomly stirred (non-rotating) turbulence [8,10]. Regarding 
the random force spectrum and statistics, we take a Gaussian random force that is white in time, for simplicity, but 
we allow for (translation invariant) spatial correlations. So we can write 



(fi(x, t)) = and (fi(x, t)fj(af, t')) = D^(x - x')8(t - t') 



(2) 



where the angular brackets denote an average over the random force realizations. The spectral function for Dij(x) 
will be specified below. 

We assume that the fluid is incompressible so that the density field is constant (p(x, t) = po) and V • u = 0. Under 
this condition we need only consider the equation for the conservation of momentum (lb) and write 



+ ^V 2 u - 2fi A u + f . 



(3) 



Notice that the force per unit mass f will be taken solenoidal as well, that is V • f = 0, in order to avoid having a 
random component in the pressure. 

In the absence of random stirring particular solutions of Eq. (3) are well known (for an incompressible fluid): it 
admits plane wave solutions, called inertial waves [1,12,13]. These are exact solutions of the non-linear equations, but 
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superposition docs not hold. They may have a role in the transition to turbulence [4]. In Subsection IV B we will 
study how the perturbative corrections modify the propagation of inertial waves. 

We now proceed to eliminate the gradient term of Eq. (3) by making use of the incompressibility condition [10]. 
We define the generalized pressure as p* Ep-^(fiA x) 2 . By taking the divergence of the previous equation we can 
solve for p* to obtain 

P* = -P0^2 ( U 3 d 3 U i) + ZtijkVjdiUk} , (4) 



so that the pressure p* can be eliminated from Eq. (3) by writing 

(u-V)u + 2ftAu . (5) 



Po V 2 



We can write the Navier-Stokes equation as follows: 

+ A7?[(u • V)u] = iA7 2 u--p(2fi Au) + f , (6) 

where, following standard practice, we have introduced the constant A in front of the advective term for book-keeping 
purposes [8] (A will be useful when carrying out the perturbation expansion and is to be set to one afterwards). The 
projection operator V is given by 

V = 1-V±V-, (7) 

and ensures that the non-linear and Coriolis terms are solenoidal. In Eq. (6), if u is solcnoidal so is f and vice versa. 

Unlike Eqs. (lb) or (3), Eq. (6) is translation invariant. That is, the centrifugal term in (lb) clearly distinguishes 
the origin (x = 0) as a special point; but as we have seen in Eq. (5) we can include this term into the generalized 
pressure and eliminate p* from the equation. This yields Eq. (6) in which a preferred direction (but no preferred 
point) is singled out by the angular velocity. This latter equation is invariant under translations, hence, we can make 
use of the Fourier transform. Since in Fourier space 

7^.(k)=Jy-^, (8) 

Eq. (6) only contains vectors orthogonal to k and we may refer to this equation as the transverse Navier-Stokes 
equation. 

We choose the random force spectrum (2) as follows: 

(fi(k,w))=0 and (fi(k,Gj)fj(k',w')) = (2D)k- y (2ir) d+1 P i:j (\i) 5(oj + lu') S d (k + k') , (9) 

where fc=|fc|,Z?>0isa measure of the amplitude of the random force and the real exponent y > — 2 characterizes 
the random force spectrum [9]. When y = d the velocity correlations produce an energy spectrum resembling the 
Kolmogorov spectrum. The random force acts isotropically on the fluid, thus, whatever anisotropies emerge at large 
scale must be due to the Coriolis force. In the following Sections we analyse the nature of the anisotropies by solving 
Eq. (6) in perturbation theory. 



B. The Fourier transformed equation 

In what follows we use the Fourier transformed equation of motion (6). Our convention for the Fourier transform 
is given by 

Uj (x,t) = f r ^ Uj (k,u,)e^-^ . (10) 

Jk<A (27T) a J-oo 27T 

We have introduced a wave- number cut-off A, so that the integral over k is restricted to the values \k\ < A. The inverse 
of this cut-off, 1/A, can be associated with the dissipation (Kolmogorov) scale, and we assume that 1/A -C L [9]. 

In order to write the Navier-Stokes Eq. (6) in wave-number representation we transform the u field according 
to (10), apply convolution to the non-linear term, and invert the Fourier transform to obtain 
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if d d p f°° dui' -> 
{-iu + vk 2 )ui(k,u) +7 3 y(k)(2f2 A u)j = --\[P lk {k)k j + 'P ij (k)k k } / — — j / — Uj(k - p, lu - iu')u k (p, lu') 

2 J P <A J-oo ^ 

+ fi(k,uj). (11) 

This equation, but without the Coriolis term, is a familiar expression in turbulence research [8-10]. Eq. (11) can be 
iterated to any desired order in A and will serve as the starting point for constructing the perturbation expansion, 
which is considered in the following Subsection. 



C. The linear response function 

If we "shut-off" the non-linear terms (proportional to A) in (11) we can identify the (inverse) linear response 
function, which will be used in carrying out the perturbative calculation. Thus, setting A = in Eq. (11) we have 

(-iu + vk 2 )u\°\k,u) +V ij (k)(2Sl Au^j =fi(k,u) , (12) 

where u' ^ is the linear velocity field. The previous equation can be written as 

[Go(fc,w)]-V°)(£,w) = f(£,w) , (13) 

which defines the linear (inverse) response function matrix, and allows one to solve for the linear velocity field 
u(°)(fc, lu). In the next Subsection we calculate the perturbative corrections to the linear inverse response function 
due to the existence of non-linearities and random forcing, and establish their effect on the physical parameters (e.g., 
viscosity) describing the transverse Navier-Stokes equation. 



D. The non-linear response function and its first order perturbative expansion 

In this Subsection we proceed to expand Eq. (11) by iteration, in powers of the non-linear coupling (see the 
diagrammatic representation, Fig. 1). In this way one can calculate the perturbative corrections to the response 
function to any desired order in A. We do not present all the details of such expansion as these can be found in review 
articles [14] and textbooks [10]. 

The calculation of the non- linear response function, [G], requires the correction matrix, [M], which is defined by a 
recursion formula: 

[G] = [G ] + [G ] [M] [G] = [G ] + [G ] [M] [G ] + [G ] [M] [G ] [M] [G ] + • • • , (14) 

which implies 

[G]" 1 = [Go]" 1 - [M] . (15) 

In Appendix A we present the diagrammatic representation of the exact transverse Navier-Stokes equation, as well 
as the first order correction to the response function. The diagram of Eq. (14) with the correction [M], calculated to 
first order of perturbation theory in the random force amplitude, is presented in Fig. 2. When written in components, 
the correction matrix [M] is given by 

M mn (fc, lu) = x 4 x [krP mj (k) + kjVmrik)] I rjn {k, lo) , (16) 

where the factor 4 is combinatorial (see Fig. 2), and the function I r j n (k,uj) is defined by 

f d d p f°° duj' 

I rjn (k,w)= ——j / — [{k s -p s )Vin{k- p) + (k n -p n )P ls (k - p)] 

Jl/L<\p\<A { Z7r > J- 

x(2D) \p\~v [G (p,cj% a V ab (p)[G (-p,-w')}l [G (k -p,cj - u%, , (17) 

(where T means the transposed matrix). The integral over p would be divergent at p — 0, so it is cut off at the low 
wave-number given by the inverse of the system size scale, 1 /L. The evaluation of Eq. (17) also involves an integration 
over the frequency a/, which has been carried out for slow rotation (see Appendix B). 
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As we are interested in the late time and large scale limits we have computed M mn {k,ui) only for u) = and up to 
second order in k [8] . The integration over wave- number p can be decomposed into an integration over modulus and 
angles. The integrand must be expanded, up to second order in k, as this is sufficient to obtain the effective viscosity 
[see Eq. (B2)]. The calculation of I r j n (k,uj) is straightforward but tedious and its technical details are presented in 
Appendix B. 

According to Eq. (15), the matrix [M] provides the order A 2 correction to the inverse response function. To obtain 
[M] (for ui = and in the limit k — ► 0), we carry out the integral of Eq. (17) and substitute into Eq. (16): 

^ A°> " -^T^WW 1 '" - y - i)k * p « (k) (18a) 

+ ^^^V-I-ZI'+V '-* 2 + " + m ' k '^ ~ , (18b) 

where Sd is the surface area of the unit sphere in d dimensions, A and 1/L are, respectively, the upper and lower 
wave-number cut-offs, A = 1, and e^-fe the Levi-Civita tensor for d — 3. In arriving at this hnal form we have assumed 
slow rotation, that is, we have kept the rotation dependent terms up to O(il) in the calculation. This is explained 
in Appendix B and we will come back to this point in Section V. We can compare the linear inverse response 
function (Bl) to its first order correction Mij(k, 0). In order to do so, we set w — in Eq. (Bl) to obtain 

[Go(k, O)]^. 1 = vk 2 % + 2e nmj Q m P m (k) , (19) 

which is split into an isotropic part (proportional to v) and an anisotropic part (proportional to 0). The first part of 
[M] (18a) corrects the molecular shear viscosity u, exactly as in the absence of rotation [8]. The non-linear inverse 
response function must take the same form as [Go(k, 0)]" 1 , with v replaced by v' , the effective viscosity [8]. The value 
of v' can be obtained by making use of Eqs. (15) and (18a), so that 



If we multiply Eq. (20) by V m i{^) and make use of the fact that V is idempotent, we obtain 

, u , D A d '^ 4 - (l/LY-y^ (d 2 -y- A)S d 

^ 2v 2 d-y-A d(d + 2)(2n) d ' 1 ' 

which shows that the isotropic term of [M], Eq. (18a), renormalizes the molecular shear viscosity v. 

We must distinguish three cases, namely, y < d — 4, y = d — 4, and y > d — 4. If y = d — 4, naive perturbation 
theory fails. Given that AL ^> 1, in the case y < d — 4, the term A d ~ s ~ 4 dominates over (1/L) d ~ y-4 , so we could 
take the limit L — > oo (the p integral is convergent in the lower limit). In contrast, in the case y > d — 4, the term 
(1/L) d ~ y ~ 4 dominates and we instead neglect the contribution proportional to A d ~ y ~ 4 . In consequence, the actual 
expansion parameter is Dk d ~ v ~ A jv 3 or DL v+4 ~ d /v 3 , according to whether y < d — A oi y > d — 4, respectively. We 
will take d — 4 < y < d 2 — 4. The number DL y+A ~ d j V 3 can be identified with the cube of the Reynolds number, 
on account of the interpretation of L as the system size scale, and that the dissipation rate is proportional to the 
amplitude of the random force D [9] . 



In particular, in d = 3, 



D_ (l/Lyy- 1 (y - 5)S 3 _ / 1 5-y DL y+1 
2v 2 y + 1 15(2tt) 3 V \ 60tt 2 y + 1 u 3 



(22) 



The most interesting case is y = d = 3, which yields the Kolmogorov energy spectrum. 

The second part of the correction (18b) has the same structure as the anisotropic term on the right-hand side of 
Eq. (19), on considering that the first term within the square brackets vanishes when it acts on Uj. However, it does 
not correspond to a renormalization of 17, owing to the additional dependence on k, namely, the k 2 factor. In fact, 
this factor is adequate for a "viscosity term" that depends on Jl (anisotropy) . Therefore, it suggests the introduction 
of an anisotropic "viscosity", which we address in the following Section. 
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III. SYMMETRY REQUIREMENTS: THE EFFECTIVE "VISCOSITY TENSOR" 



In the previous Section we have seen that the term (18b) induces no corrections to v. In this Section we give a 
physical interpretation of this anisotropic contribution, and show that it arises from an effective "viscosity tensor" for 
the rotating fluid with mixed symmetry in its indices. For a Newtonian fluid the linear relation between the rate of 
strain and stress tensors involves a rank four viscosity tensor, so that we can write [15] 

_ 1 / du m dun \ = . . 

1 ij ~ r,Vijmn r. + a I — VijmnU mn ■ 



2 lljmn \dx n ' dx n 

As both Tij and u mn are symmetric tensors the only symmetries of the "viscosity tensor" are the following: Vijmn — 
Vjimn and rjij mn — Vijnm- From the above symmetries we conclude that the "viscosity tensor" has 36 independent 
components (in d = 3). We write Vijmn as a sum of a pair-symmetric (S) and a pair-antisymmetric (A) part as follows 

Vijmn = ^{Vijmn + Vmnij) + ^(Vijmn — Vmnij) = Vijmn + Vijmn ) (24) 

so that vfjmn has the same symmetries as Vijmn plus Vijmn = iLij > an d Vijmn nas ^ e same symmetries as Vijmn phis 
vtjmn = ~Vmnij- There are 21 pair-symmetric and 15 pair-antisymmetric independent components. The presence of 
axial symmetry (induced by the Coriolis term) reduces the number of independent components of both vfjmn an( ^ 
vtjmn- The most general axisymmetric tensor (in d = 3) can be constructed from Qi, 5ij, and e^fe as follows 

Vijmn = a 1 {fl 2 )(S tm S jn + SinSjm) 

+ a 3 (Q 2 )(H, i Q : jS mn + flm^nSij) 

-\- 0^4(^1 ) (ilii^lm^jn ~t~ ^j^m^in ^i^n^jm ~t~ ^j^n^im) 

+ a5(£l 2 )Sljf2j£l m f2 n , (25a) 



— )^q{^qim^jn ~t~ £qin&jm ~\~ ^qjm&in ~t~ ^qjn^im) 

~\~ ^2(^2 )^q^qim^j^n ~1~ ^qin^j^m ~t~ ^qjm^i^n ~t~ ^gjn^i^m) 



+ 3 (fi 2 )(fi j fi J -j inn - fi m r>„%) . (25b) 

We observe that the number of independent components has been reduced from 21 to 5 for the pair-symmetric term 
and from 15 to 3 for the pair-antisymmetric part. The symmetry arguments used in deducing (25a, 25b) are analogous 
to those used in the theory of elasticity [16]. 

The coefficient functions ai(Q 2 ) and Pi(£l 2 ) can be written as a series in f2 2 : 



+00 +00 



n„ = % a nr (n 2 y , p n = J2Pnr(n 2 ) r . (26) 



In the limit of slow rotation (linear order in f2) these coefficients reduce to their constant (il = 0) value and, 
furthermore, the terms proportional to 0:3, 04, a 5 , /3 2 , and /?3 vanish at this order. We can then write for the "viscosity 
tensor" 



Vijmn = ai(SimSj n + 5 in S jm ) + a 2 &i 3 5 mn , (27a) 

Vijmn ~ Pl^qi^qim^jn ^qin^jm H~ ^qjm^in ~\~ ^qjn&im) • (27b) 

The pair-symmetric part is isotropic and we can identify poK = a 2 + 2a\/d and p$v — cti, where k and v denote 
the molecular bulk and shear kinematic viscosities, respectively; the pair-antisymmetric part, proportional to the 
coefficient is identified below. Incompressibility means that n — 0, leaving only v, in the absence of rotation. 
This is the molecular viscosity in the original Navier-Stokes equation (lb), which gets renormalized, becoming an 
effective (rotation independent) viscosity, as seen in the previous Subsection II D. However, it is to be expected that 
perturbation theory, at sufficiently high order, generates a dependence of ct\ on ft, such that one would be led to 
consider a rotation dependent viscosity. At the order we are working, the effect of rotation is to generate a rotation 
dependent anisotropic "viscosity" of the type (27b), as will be demonstrated below. 

The anisotropic "viscosity" (27b) is associated to a stress tensor T£, which yields the following viscous force 
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F t = d 3 T* = d 3 ( V f 3mn u mn ) . (28) 
In wave-number representation this can be written as 



Fi — ^-^ kjr/^ mn (k n S mq + k m 5 nq ) u q — kjk m rj. 



A 

ijqm ' 



u q , (29) 



where we have made use of the symmetry within pairs of indices of the tensor r\ A . In order to be consistent with the 
transverse Navier-Stokes equation, we must project Eq. (29) by means of 'P(k) 

V v (k) F p = -kjkmVMv^ u q . (30) 

We now substitute Eq. (27b) into the previous expression to obtain 

^zp(k) Fp — /?i7^p(k) kj kniQni^tpqnfijm £pmn&jq £jqn$pm ~t~ £jmn$pq) Uq 
— Pi \kqkm£imn 

+ k l ^mgn^im(k)] U q — {3\Q<nk tmqn'Pimfit) Uq : (^-Q 

where we have taken into account the antisymmetric properties of the Levi-Civita tensor and the fact that Vij(\i)ki = 
Vijitykj = 0. We now compare this equation with the first order perturbative correction of the linear inverse response 
function. The force given by Eq. (31) agrees identically with the force provided by the anisotropic part of the correction 
Mij(k,0), namely, the product of its expression in (18b) and Uj, once we choose the coefficient in (27b) to be 

_ (2D)(2S d ) 2 AW-(1/L)W 
Pl = - P0 ( 27 r)W ( ] (d-y-6)d(d + 2) " (32) 

We have thus shown that the contribution (18b) to the response function can be understood as arising from the 
anisotropic "viscosity" tensor (27b). Since we are considering the case y > d — 4, we can neglect the term A d ~ v ~ 6 . 
The expansion parameter, proportional to D/v 3 , must again be identified with the cube of the Reynolds number. Wc 
postpone this identification to the following Section. 
If we set d = 3 in the expression for /?i , wc obtain 

(2D)(2S 3 ) (l/L)-v-z p DLV+ 3 

Pl P0 (27r)3(2^)3 (y + 3)15 157r 2 (y + 3) ^ ' [ ' 



IV. EFFECTIVE LARGE SCALE DYNAMICS 

In this Section we focus on the new force Fi — djT^ that appears in the fluid equation of motion due to first order 
perturbative corrections, and discuss some possible physical implications. 

A. The quasi-local force 

Here we use the results derived in previous Sections to "correct" the transverse Navier-Stokes equation with the 
newly generated (first order in fi) terms that arise at large scales due to rotation. In coordinate space the force (28) 
becomes 



Fi=0i 



(n A V 2 w) 4 + di[Q ■ (V A u)} . (34) 



In deriving this equation we have made use of the incompressibility condition and the antisymmetry properties of 
the Levi-Civita tensor. This is a quasi-local 1 force (per unit volume) in which the angular velocity couples to the 
Laplacian of the fluid velocity and to the vorticity u = VAu. 



1 Quasi- local, in general, means depending on a function (in this case the velocity field u) and its derivatives. 
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The physical character of the force Fi can be revealed by writing the effective equation of motion in coordinate 
space. We have 

^ + XV(u • Vu) = i/V 2 u - V(2Sl An) — P[ V(2Sl A V 2 u) + f , (35) 

where v' is the rotation independent effective kinematic viscosity, [see Eq. (21)], and f3[ — (3i/(2p ). The magnitude 
of the correction linear in f2 (the new force F) can be determined by comparing it with either the Coriolis force or 
the viscosity correction. We obtain for the ratio to the Coriolis force: /3[ |(V 2 u)| / |u| ~ DL y+1 /v 3 (LA) 2 , assuming 
that the scale of spatial velocity fluctuations is 1/A, that is, |(V 2 u)| / |u| ~ A 2 . Hence, the expansion parameter can 
be identified with the cube of the Reynolds number (as in the previous Section) times (LA) 2 . On the other hand, the 
ratio between the correction linear in fi and the correction linear in v' is ~ VlL 2 jv ~ Ek^ 1 . This was to be expected, 
since Ek measures the relative strength of the viscosity and Coriolis force. 

We conclude this Subsection by providing an important property of the new force. In general, for a stress tensor 
associated with a pair-antisymmetric "viscosity tensor", the power is given by [13] 

P oc J d 3 x UijTfi = J d 3 x Uijr]^ mn u mn , (36) 

but, as f]ij mn = —Vmnij ! we conclude that P — 0; therefore, this stress tensor, T-j, does not lead to dissipation and 
is not truly viscous. This implies that the name "viscosity tensor" is not appropriate, and we have only introduced it 
by analogy with the truly viscous pair-symmetric tensor [13]. 



B. Inert ial waves 



As already mentioned, rotating incompressible fluids support wave solutions that are exact solutions of the non- 
linear equations [1,12]. Although we have derived the quasi-local force from the randomly forced equations (that 
describe turbulence), it is of interest to see how this force affects these wave solutions. Let us consider wave solutions 
of the form e I ( fc :E + wt ) as single mode plane waves of the linearized equation (35) (the incompressibility condition 
annihilates the advective term) but without the forcing term 

^ = t/V 2 u - V(2n A u) - /3[P(2n A V 2 u) . (37) 

We can write this equation in components as follows 

(iw + v'k 2 ) Ul = 2(-l + /3ifc 2 )P„„(k)e m9 „f> 9 M„ . (38) 

We are free to choose f2 = (fl x , 0, £l z ) and k = (0, 0, k). This choice implies that u z = (incompressibility condition), 
so that the waves are transverse. From our choice for k the only non- vanishing components of 'P(k) are V xx and V yyi 
which are equal to one. We can write the x and y components [the z component of (38) is identically null)] as follows 

(iuj + v'k 2 )u x = -2(-l + x k 2 )Sl 7 u y , (39a) 
(iuj + v'k 2 )u y = 2(-l + l3[k 2 )Q z u x . (39b) 

These equations yield the (complex) frequency of the plane waves (dispersion relation) 

u(k) = ±21} cos (9(1 -f3[k 2 ) +iv'k 2 , (40) 

where 9 is the angle between k and fl, so that k • fl = kQcosO, with k — \k\ and ft — |f2|. By making use of the 
dispersion relation (40) it is easy to see that the waves are circularly polarized and transverse 

u y = ±iu x , u z = . (41) 

Wave packets are not solutions of the non-linear equations and can only be considered for small wave amplitudes. 
The group velocity V of a wave packet is (for vanishing viscosity) 

VS) = ^ = ± ^Zm vMaj T 4# £n • k , (42) 
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and the phase velocity v (for vanishing viscosity) 



k t = ±2(1 - (3[k 2 ) 



ft k 



(43) 



The "standard textbook" result is recovered by taking the limit (3[ ^ [1,12]. From this calculation we see that the 
new term in Eq. (37) not only changes the wave frequency and phase velocity, but also the group velocity. Moreover, 
the group velocity is no longer perpendicular to the phase velocity 



The quasi-local force would cause the energy transport to not be perpendicular to the phase velocity and a small 
fraction of the energy in the waves to be transported parallel to the wave vector. 



We have applied perturbation theory to a homogeneous incompressible viscous fluid subject to solid body rotation 
and isotropic random forcing. At small scales, we assume that only the molecular shear viscosity and the Coriolis 
force are needed in writing down the transverse Navier-Stokes equation (6). Our (first order) perturbative results 
demonstrate that (i) anisotropic components of the effective "viscosity tensor" are dynamically generated at large 
scales by the combined interplay of the Coriolis force, the random forcing term, and the inherent non-linearity of 
the Navier-Stokes equation, (ii) the molecular shear viscosity v gets corrected in the same manner as for isotropic 
randomly stirred turbulence [8-10]. 

These perturbative results are corroborated by a symmetry principle. By making use of the axial symmetry of 
a rotating fluid, we have constructed the most general "viscosity tensor" that is invariant under such symmetry. 
The preferred direction singled out by the angular velocity ft breaks the isotropy and leads to new terms in the 
"viscosity tensor" absent in the isotropic case [13]. We have also determined and described the role of the first order 
perturbative new term in the effective fluid equation of motion. It acts as a quasi-local force and, like the Coriolis 
force, is not dissipative and does no net work on the fluid. Most importantly, we find that this quasi-local force 
affects the propagation of inertial waves in rotating fluids. For small wave amplitudes a fraction of the wave energy 
is transported in the same direction as the phase velocity. 

The perturbative calculation is developed as a double expansion: in both the amplitude of the random force and 
the angular velocity. The actual dimensionless expansion parameters turn out to be the cube of the Reynolds number 
and this number times the inverse of the Ekman number, respectively. For simplicity, we have restricted ourselves to 
the computation of the lowest order in both: that is, first order in the random force amplitude and linear order in ft. 
This is sufficient to generate two, out of the total of eight tensor terms (seven, on account of incompressibility) allowed 
by axial symmetry [see Eqs. (25a) and (25b)]. We conjecture that, at first order in the random force amplitude, all 
the remaining tensor terms are generated for higher powers in fl: 013, and f3% at quadratic order, P2 at cubic 
order, and finally a$ at quartic order. Of course, as increasing powers of fl are taken into account, the coefficient 
functions (26) must be expanded out to the order of ft being investigated. The higher order terms will be needed to 
study the effects of fast rotation and to track the onset of the bi-dimensionalization of the fluid [4] . 

An important step in this direction will be provided by a complete renormalization group (RG) analysis of the 
large scales properties of a rotating incompressible fluid. In order to carry out this RG analysis one may need to 
calculate the (perturbative) corrections to the non-linear coupling term A, as well as the random force amplitude D, 
and combine these two with the response function calculation presented in this paper. Once in hand, the RG fixed 
points can be determined and the corresponding asymptotic behavior of the rotating fluid deduced, allowing us to 
compute quantities such as the scale-dependent Reynolds and Ekman numbers, among others. We hope to report on 
these developments elsewhere. 



The authors thank Juan Perez-Mercader and Alvaro Dommguez for discussions and for reading an earlier version 
of this paper. One of us (DH) acknowledges correspondence with Arjun Berera covering a wide range of topics in 
turbulence theory. Our work is supported by grant BFM2002-01014 and the work of Jose Gaite is further supported 
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V-k = T4/3ifcr2-k^O. 



(44) 



V. DISCUSSION 
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APPENDIX A: DIAGRAMMATIC REPRESENTATION OF THE EXACT EQUATION AND THE FIRST 

ORDER CORRECTION TO THE RESPONSE FUNCTION 



In this Section we include a diagrammatic representation of the exact transverse Navier-Stokes equation [see Eq. (11)] 
and of the recursion relation for the first order response function [see Eq. (14)]. We have followed the convention 
introduced in Ref. [8]. 

The exact equation can be represented as 




FIG. 1. Exact hydrodynamic equation. 



where the thick lines represent u(k, u>), the open circles represent f(k,w), and the thin lines represent Go(k,w). This 
graphic representation implies that the linear velocity field u^°\k, u>) is given by a thin line attached to an open circle, 
as we have u( ) (k,w) = Go(k,u)f{k,u). 

In this diagrammatic representation we can also depict the recursion relation for the response function Gij(k,u), 
Eq. (14), given by the first order correction [M], [Eq. (16)], as follows 




k- p 

FIG. 2. Recursion relation for the first order response function. 

where the solid circle represents the random force average (f m (p, w')/ n (— P, — w')) introduced in Eq. (9), and the 
double thin lines represent the first order response function [G]. 

APPENDIX B: PERTURBATION EXPANSION AND EVALUATION OF THE INTEGRALS 

The large distance and long time renormalizability of the transverse Navier-Stokes equation (6) implies that the 
corrected response function G must have the same mathematical structure as its linear counterpart Go- G has, 
therefore, the same frequency and wave-number dependence as Go and contains the same number of parameters. The 
renormalization of the response function yields the (rotation independent) effective viscosity v' . In this Appendix we 
outline the major points in the calculation of the first order correction to the response function [8]. 

Eq. (13) defines the linear inverse response function, which is given by 

fa-b h 0\ 

[Go(fc^)]- 1 = / a + b , (Bl) 
\ c da) 

where the matrix entries are: a = —iu + vk 2 , b = 2flkik 2 /k 2 , c = —2Qk 2 k 3 /k 2 , d — 2£!fc 1 fc 3 //c 2 , / = 20(1 — k%/k 2 ), 
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and h = — 20(1 — kf/k 2 ). In the limit of vanishing rotation (O = 0) we recover an isotropic diagonal matrix for the 
linear response function and its inverse [8] 

[G (k, lu)}^ 1 = (-iuj + vk 2 ) % , [G (k, = (-iw + vk 2 Y x 5 t] . (B2) 

Given the matrix form of the inverse linear response function (Bl) we can write for the linear response function 2 

, / a 2 +ab -ah \ 

[G (k,cj)} = - 3 -af a 2 -ab . (B3) 

a ab a J h \-ac-bc + df -ad+bd + ch a 2 -b 2 -fhj 

The first order correction to G is given in Subsection II D [see Eqs. (14), (15), (16), and (17)]. We have restricted 
ourselves to the limit of slow rotation (linear order in f2), and therefore need to expand both Go and Gq 1 to this 
order. We point out that Gq 1 is already linear in ft [see Eq. (Bl)]. To linear order in £1 the linear response function 
is given by 

1 /a + b -h 0\ 1 /l 0\ , / b -h 0\ 

[G (fc,w)] = -2 -/ a- b +0(fi 2 ) = - 1 +— -/ -b 0)+O(n 2 ). (B4) 
a V -c -d a) ° \0 l) a V-c -d 0/ 

The first step of this calculation is to carry out the frequency integration (over a/) in Eq. (17). From Eq. (19) we see 
that the viscosity v is the coefficient of the k 2 term, and to obtain the effective viscosity we can set the frequency 
ijj to zero from the outset. Once we make u> = the integral over oj' can be computed by means of the calculus of 
residues [8]. We close the contour in the lower half plane, where there are in general three simple poles, which coalesce 
into one double pole in the limit of slow rotation. We obtain (keeping up to linear terms in ft): 



f 

j — ( 



du)' (1 + ^r) (1 + 2i? 2 fc ) 

— [G (p, u')]rt[G (-p, -U% q [G (k - p, = S rt 5 sq Sjl ^^2)2 + Srt Ss 9 {2 V p^f 



(1 _|_ £^) (2+2^) 

+ 5 rt Sji (2 ^ 2 p j 3 M-flU + tsgSfl (2 ^ 2 p j 3 [n(p)} rt 

+ 0(k 2 ) . (B5) 

We have carried out an expansion in k = \k\ and only kept up to linear order, as Eq. (16) is already linear in k, and we 
only need to compute the k 2 term. The previous equation was evaluated assuming a positive molecular shear viscosity 
coefficient {y > 0). A change in sign will change the location of the poles and modify the frequency integration. We 
have introduced the matrix [fi(p)] defined as follows 

[Q(p)]mk = -2e ijk 9.j V im (p) , (B6) 

which is already linear in the angular velocity. 

If we make use of this intermediate result (B5) and substitute it into Eq. (17), we are left with the integration 
over the wave-number p. From Eqs. (15) and (19) we can see that the renormalization of the molecular viscosity v 
requires that we expand the first order correction [M] up to second order in the wave-number k. The factor in square 
brackets in Eq. (16) is already linear in k, so that the integral (17) only needs to be expanded to first order in k. 
It is important to note that this integral depends on k not only through the integrand but also through its limits 
of integration. This is because all wave- numbers appearing in (17) must remain within the set of wave- numbers to 
be integrated over (in this case we must ensure that both p and k — p belong to this set) [8] . This means we must 
integrate over the intersection of the domains 1/L < \p\ < A and l/L < \k — p| < A, where 1/L is the lower cut-off. 
To first order in k the second inequality can be written as 1/L + k cos 9 < p < A + k cos 9, where 9 is the angle between 
k and p (k ■ p — kp cos 9). There are two cases to consider: (i) if cos 9 > the intersection of the two intervals can 
be expressed as the difference of intervals [l/L, A] — [1/L, 1/L + kcos9] and (ii) if cos 9 < the intersection can be 
written as [1/L, A] — [A + kcos9, A]. This means that the complete wave-number integration, valid up to 0(k 2 ), can 
be written as 



2 This is also sometimes called the linear propagator [8,10]. 
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where d£l d is the surface element of the unit sphere in d dimensions. 

We also present the following projection operator product expansions that prove to be useful when handling the 
intermediate steps of the computation 

v aM {v V )v dc {v) = (K -p a )(5 bd - M*) _ E^(k d _ tip.} + (k 2 ) , 

7> Q , bc (k - p)T ad (p) = k a T bc (p)V ad (p) + Oik 2 ) , (B8) 

where 

V a , bc {k) = k b Vac(k) + k c V ab (k) . (B9) 

Finally, various identities needed for the angular integrations are collected here. Let S d represent the surface area of 
the unit ci-sphere and hj denote a unit vector in the j-th direction. The only angular integrations required are of the 
following types 

J dfl d 
J dfl d hihj 
J dfl d hihjh n h m 

The angular integration of the product of an odd number of unit vectors over the unit d-sphere vanishes identically. 
If we make use of the previous results, we obtain the first order correction as written in (18a, 18b). 



S d , 
d ° 13 ' 

did + 2) ( ( ^^ mn ^ im ^3 n + SinSjm) . (BIO) 
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